The linear sigma model at a finite isospin chemical potential 
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The effect of finite isospin chemical potential to the effective masses of the mesons at finite 
temperature is investigated in the framework of the 0(4) linear sigma model with explicit chiral 
symmetry breaking. We present a mechanism to include the isospin chemical potential in the 
model. By using the Cornwall- Jackiw-Tomboulis method of composite operators, we obtain a set 
gap equations for the effective masses of the mesons and get the numerical results in the Hartree 
approximation. We find that the introduction of the chemical potential only affects the mass of the 
charged pions and sigma, while there is almost NO effects on the mass of neutral pions. 
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I. INTRODUCTION 



Nowadays, it is widely believed that QCD is the most promising theory to explain the strong interactions. However, 
due the property of confinement and asymptotic freedom we can accept and study the quarks as free only at very high 
energies. This much expected phase, which is called Quark Gluon Plasma (QGP), is under investigation in RHIC 
experiments and there is strong evidence that it has been observed. According to Heinz pj it seems that elliptic 
flow at high enough collision energies is a QGP signature. The recent results from RHIC experiments have attracted 
many attempts for theoretical models to be testedQ- Recently there are also a lot of attempts to study the phases of 
\Q . QCD from first principles using lattice techniques at finite baryon density and chemical isospin potential. There is, 
however, an alternative approach which applies only in the low energy regime. Following this path, it is necessary to 
use a variety of effective phenomenological models and QCD-like theories in order to mimic some of the features of 

C^; QCD. 

A quite popular model, which fulfils these requirements, exhibiting many of the symmetries observed in QCD, is the 
J^" 1, linear sigma model. The model was first introduced in the 1960s as a model for pion-nucleon interactions Q and has 
attracted much attention recently, especially in studies involving disoriented chiral condensates 0, 0, 01 ■ The linear 



sigma model with a U(Nf)x U(Nf ) symmetry and two to four quark flavors have been studied in the Hartree- Fock 



or Hartree approximation |7l l8l lii ll(| within the Cornwall-Jackiw-Tomboulis (CJT) formalism [llL and the large-N 
approximation of the linear sigma model has been investigated by using the same formlism in Refs|l2l Within 
' the same framework, Shu and Li |15||have studied Bose-Einstein condensation and the chiral phase transition in the 
chiral limit, and they introduced the chemical potential through the third component of isospin charge. In this 
discussion, we apply a mechanism [Lsl Il7l Il8j | to include the isospin chemical potential in the linear sigma model, 
then by using CJT method, we study the chiral phase transition at finite isospin chemical potential in explicit chiral 
symmetry breaking. 

The organization of this paper is as follows. In the next section we present the basics about the linear sigma model 
and how this model incorporates some of the basic features of QCD, then we introduce the isospin chemical potential 
in the model. In Section HTTl we calculate the effective potential of the model by using the Cornwall-Jackiw-Tomboulis 
method of composite operators and obtain a set of gap equations for the thermal effective mass of the pions. The 
solution of these gap equations is presented in Section llVI In the last section we present our conclusions. 
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II. THE MODEL 



Consider an idealisation of QCD with two species of massless quarks u and d. The Lagrangian of strong interaction 
physics is invariant under SU{2)l x SU(2) r chiral transformation 

*z,fl -> ex.p(-ie L , R ■ t)^l,r, (1) 

where ^l.r — {u,cI)l.r- However this chiral symmetry does not appear in the low energy particle spectrum since 
it is spontaneously broken. Consequently, three Goldstone bosons, the pions, appear and the (constituent) quarks 
become massive. At low energy, the spontaneous breaking of chiral symmetry can be described by an effective theory, 
the linear sigma model, which involves the massless pions tt and a massive a particle. Here, the a field can be used 
to represent the quark condensate which is considered as the order parameter of the chiral phase transition. This 
choice is a natural one since both exhibit the same behaviour under chiral transformations. The pions are very light 
particles and can be considered approximately as massless Goldstone bosons. 

There are different ways to write down the Lagrangian of the linear sigma model and each with certain advantages 
and disadvantages. One version which demonstrates the chiral properties of this Lagrangian in a straightforward way 
involves the introduction of a matrix field $ defined as 

T° f 

$ = (7 _ +Z7? ._ j (2) 

where r° is the unity matrix and f — (ti,T2,Ts) are the Pauli matrix. This matrix field satisfies the normalization 
condition Tr(r a T ) — 2S ab . Under SU(2)l x SU(2) r chiral transformations, $ transforms as 

$ -> L+$R . (3) 

Using this parametrisation the renormalizable effective Lagrangian of the linear sigma model can be written as fl9| 

m 

£o = U + C q , (4) 

where 

= Tr[(<9 M $) + (<^$)] - m 2 Tr($+$) - ^[Tr($+$)] 2 (5) 



C q = * L vfdyp L + ^Rvfd^R - 2gV L $* R + h.c. (6) 

When chiral symmetry is breaking, the field a acquires a non-vanishing vacuum expectation value, which breaks 
SU(2) L x SU{2) R down to SU{2) L+R . It results in a massive sigma particle a and three massless Goldstone bosons 
7r, as well as giving a mass m q = gf n to the constituent quarks. 

Now let us consider the Lagrangian of the linear sigma model for two massless quarks of flavors u, d with different 
chemical potential fj, u and fi c i, 

C' = £o+£» = U +C q + y{BV L + tf^B* R) (7) 



where B = 




is the matrix of chemical potentials. The term £ M can be expressed either by the variables 



M«) or by the two combinations fi R — hif^a + Md) an d Ml — i(M« ~ ^d)> which couple to the baryon charge density 
and the third component of isospin respectively. Then we have 

£ M = ,i B (*[* L + tf R nr R ) + + • (8) 

Both and C q are SU(2)l X SU(2) R invariant, but the symmetry is reduced by the baryon charge density term 
fj, B (^L^ L + ^r^r) to SU(2) l +r and further reduced by the isospin term ^/(*[r 3 * L + ^ r t 3 ^ r ) to U(l) L+R . 
In what follows, we consider an idealized case where /i/ is nonzero while [i R = 0. A similar choice is adopted in 
Refs.0][2l]]. According to their analysis such an approximation is justified by the fact that we are dealing with the 
dynamics of strong interactions. In any realistic setting, /i/ <C /is- Such a system is unstable due to weak decays 
which do not conserve isospin. However, we can still imagine that all relatively slow electroweak effects are turned 
off. 
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We note that the effects of the isospin chemical potential /ij can be included to the Lagrangian of the linear sigma 
model without additional phenomenological parameters by promoting_S'?7(2)i x SU[2) R to a local gauge symmetry 
and viewing m as the zeroth component of a gauge potential 0, ITilllSj] . This means that we can rewrite the terms 



Cq + as 



£' = *Li7"A.*L + ^ R iY l D^ R , (9) 



where = — ifiiA^, — A^T a and T a = ^(a = 1, 2, 3). Here the gauge fields A^ have the following properties: 
= l,Al = Al = and A t . = for i = 1,2,3. Then the term £^ is invariant under the SU(2) L x SU{2) R local 
transformations: 

-> U9 LtR , U = exp[-i6 L>R (x) ■ f) . (10) 

However, the term is not invariant under SU(2)l x SU(2) r local transformations. In order to ensure the effective 
Lagrangian of the linear sigma model is also invariant under this local symmetry we have to replace the derivatives 
in Eq. (jSJ) by the covariant derivatives: 

U = Tt[(D^)+(D^)} - m 2 Tr($+$) - ^(Tr($+$)) 2 , (11) 
where the covariant derivatives are defined as: 

3 3 

D $ = d $-z/i/(y$-$y) (12) 
A:$ = di$,(i = 1,2,3) . (13) 

Thus the gauge invariance fixes the way fj,i enters the mesonic sector of the effective Lagrangian of the linear sigma 
model. 

In order to study the effects of the finite chemical potential fij in the linear sigma model at finite temperature, we 
only consider the term neglecting the fermion fields. It is convenient to define the new fields: 

, 7T 1 ± Z7T 2 , , 

* - ^pr - (14) 

The Lagrangian £$ now can be expressed as 

C = l(d,*)(d»<r) + ±(d^°)(d»ir°) + (D^)+(D^)- - V[<r,n) , (15) 



where the potential is 



V(a,n)^^+^) + ^ + nr , (16) 



and (D tt) ± = [do ± m)^, {D l n) ± = 5 1 tt ± for [i = 1, 2, 3). 

The SU[2)l x SU(2) r symmetry of the linear sigma model is explicitly broken if the potential V(a, 7?) is made 
slightly asymmetric by adding to the basic Lagrangian in Eq. Q a term of form 

C b = ea . (17) 

With this addition, the vector isospin SU(2) symmetry remains exact, but the axial SU(2) transformation is no longer 
invariant. The potential now has the form 

V {a, t?) = ^(a 2 + tt 2 ) + A (a 2 + -*)2 _ ea (18) 



2 v ' 24 

)1 tA „ - f J 



To the leading order in e, this shifts the minimum of the potential to v = f v + -Mr, an d as a result the pions acquire 



At tree level and zero temperature the parameters of the Lagrangian are fixed in a way that these masses agree with 
the observed value of pion mass m v = 138 MeV and the most commonly accepted value for sigma mass m CT = 600 
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FIG. 1: Two-particle irreducible graphs which contribute to the effective potential of a \(j) 4 theory in the CJT method up to 
the three-loop level. The solid line represents the dressed propagator G(x,y). There are two kinds of vertices: a three-particle 
vertex and a four-particle vertex. 

MeV. Then, the coupling constant A of the model can be related to zero temperature properties of the pions and 
sigma through the expression 

X= 3(ml - ml) , (19) 

J 7T 

where /„• w 93 MeV is the pion decay constant. The negative mass parameter m 2 is introduced in order to ensure 
spontaneous breaking of symmetry and its value is chosen to be 

-m 2 = (m 2 - 3m 2 )/2 > . (20) 



III. THE CJT EFFECTIVE POTENTIAL 



In this section we calculate the finite temperature effective potential for the linear sigma model at finite isospin 
chemical potential. In our calculations we use the imaginary time formalism, which is also known as the Matsubara 
formalism |22l I23L I24I I25J ] . This means that in the case of Bosons, we have 

d 4 k lr^f d 3 k 



where /3 is the inverse temperature, {3 = jr^T' anc ^ as usua ^ Boltzmann's constant is taken as ks — 1, and uj n = 2TrnT, 
n = 0, ±1, ±2, ±3, . . . For simplicity we have introduced a subscript [3 to denote integration and summation over the 
Matsubara frequency sums. 

The starting point is the Lagrangian given in Eq. I|15|l with the choice of the asymmetric potential given in Eq. JTSJ. 
By shifting the sigma field as a — > a + <j> the classical potential takes the form 

tf(0 = i m V + ^ 4 -^. (21) 

The tree-level sigma and pion propagators corresponding to the above Lagrangian have the form 

D- 1 = ojI + k 2 + m 2 + ^\(j> 2 , (22a) 

D^ 1 = ojI + k 2 + m 2 + ^\6 2 , (22b) 

D' 1 = (u n + i/x/) 2 + k 2 + m 2 + \\6 2 . (22c) 

D 

Eq. (|22c() represents the inverse propagator of 71+ and 7r_ due to the facts that uj n + and w„ — ifij are equivalent 
in describing the propagators of n+ and 7r_ ■ 

The interaction Lagrangian which describes the vertices of the shifted theory is given by 

r ^ 2^2 ^4 ^^4 ^ / -a 1 3 mo\ 

Li n t = <7 7T a 7T 6(771 (bcj . 23 

12 24 24 6^6^ K ' 

Here, we omit the constants and the terms linear in the a fields for simplicity. In our approximation we do not 
consider interactions which are given by the last two terms in the above Lagrangian. 
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a TT TT 

(a) (b) (c) 

FIG. 2: The Hartree contributions to the CJT effective potential. Solid lines correspond to G a , while dashed lines correspond 
to G n - 



Then, for constant 
given by 



according to CJT formalism [Tl| the finite temperature effective potential for X<p A theory is 



V{4>,G) = u{4>) + - 



lnG^O; k) 



+ 2 J [D-HfrVGi frk) 
+W,G). 



(24) 



Here, G(4>\ k) denotes the full propagator of the theory and we have introduced a shorthand notation for the space 
integration and summation over the bosonic Matsubara frequencies so the symbol f„ stands for TJ2 n J d 3 k/(2Tr) 3 . 
The last term V2(<j>, G) represents the sum of all two and higher-order loop two-particle irreducible vacuum graphs of 
the theory with vertices given by £ int and propagators set equal to G{<j>; k). The diagrams contributing to V^fy, G) for 
X<p A theory are shown in Fig.^ Evaluating the effective potential in Hartree approximation means that one needs to 
take into account the "oo" type diagram only. In the case of X<p A theory this diagram is the leading order in 1^(0, G) 
in both the loop expansion and the 1/N expansion. We adopt this approximation in the case of the linear sigma 
model as well. The Hartree and large N approximations at zero chemical potential have been studied previously in 

Refs. |1 IE II G3 [II GS [13 . 

Generalising this result to finite isospin chemical potential, we derive the following expression for the effective 
potential at finite temperature 



V{ft>,G) = U^ + ^JlnG-^^ + ^JlnG^ifrk)- 

r {Dn 1 (frk)G (frk) 

* J (3 

+ I [D-\(j>; fc)G c (</>; k) - 1] + V 2 {<p, G a , G , G c ). 



+- I [^- 1 (0;fc)G ff (^;fc) 



11+ - 



\nG-\4>;k) 
1] 



(25) 



The first term U(<p) on the right handside is the classical potential and the last term V2(</>, G a , Go, G c ) denotes the 
contribution from two-particle irreducible diagrams. In the following we include only the two-loop diagrams shown 
in Fig. These contribute the following terms in the potential 



V 2 {<l>,G a ,G v0 ,G c ) = 3 



24 



I 2 



+4 



24 



+4 



G-\cf>;k) 
G-\frk) 



X 







24 

G- l {<t>; k) 



G-\4>;k) 







24 



Go\cj>; k) 



G-\^k) 



X 
^24 



G-\^k) 



(26) 
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Minimizing the effective potential with respect to full propagators we obtain the following system of nonlinear gap 
equations: 

G- 1 = D^ + ^J G CT (</>; k) + ~J G o (0; k) + - jf G c (0; k), (27a) 

G - x = J D- 1 + ^G ff (0;fc) + ^G o (0;fc) + ^G c (0;fc), (27b) 

G- 1 = D^ + ^J G a (</>; fe) + ~ y G o (0; fc) + y jT G c (0; fc). (27c) 

In order to proceed we can make the following ansatz for the full propagators: 

G- 1 = ul + P + M 2 al 



Go 1 = L0 2 n + k 2 + Ml 
G- 1 = K + z/i/) 2 + k 2 + M t 



T-2 

C 1 



where M ai Mq and M c are the effective masses of a meson, ttq and charged pion dressed by interaction contributions 
from the diagrams of FigEl The dressed sigma and pion masses are then determined by the following gap equations: 

M 2 a = m 2 + ^ 2 + ^F{M a ) + ^F(M ) + ^F(M C ), (28a) 



2" 2 v '6 v ' 3 
^ 2 + ^F(M CT ) + ^F(M ) + ^ 



M 2 = m 2 + — 2 + —F(M a ) + —F(Mq) + —F(M C ), (28b) 



A/2 = m 2 + ^ 2 + ^(M^) + ^(Mq) + ^■F(M C ), (28c) 



6 T 6 v ' 6 v 3 
Here we have used a shorthand notation and introduced the function 



F(Ma) = I r \ — , (29a) 

J/3 w 2 + k 2 + Ml 

F(Mo) = I - * (29b) 

7/3 ^2 + fc2 + M 2 

F W = / ~, . J f2 „2 - (29c) 

J/3 (u fl + z/i/) 2 + fc 2 + Ml 

At the level of the Hartree approximation, the thermal effective masses are independent of momentum and are 
functions of the order parameter <fr and the temperature T. Moreover, in our case it is also function of the finite 
isospin chemical potential /i/. 

In terms of the solutions of the gap equations (|28|l . the effective potential at finite temperature can be written in 
the form 

VfaM) = \m 2 4> 2 + ^4> i ~^+\J\nG- 1 + \J\nGi ) 1 +J\nG- 1 

-\{Ml -m 2 - \\4, 2 )F{M a ) - 1 -{M 2 ~m 2 ~ ±\4, 2 )F(M ) 

-(Ml - m 2 - \\<j) 2 )F{M c ) + ±IF(M„)} 2 + ±{F(M )} 2 + ±{F(M C )} 2 

+±-F(M a )F(M ) + ^F(M a )F(M c ) + ^F(M )F(M C ). (30) 
12 o o 

In above expression and in following discussions, we have taken the symbols M a , Mq and M c just as the solutions of 
the gap equations (|2*5|l for simplicity. 

Minimizing the effective potential with respect to the full propagators, we have found the set of nonlinear equations 
for the effective particles masses which are given by Eq. (128(1 . In addition, by minimizing the potential with respect 
to the order parameter <fi we obtain one more equation 

[™ 2 + \^ + ^F(M a ) + ^F(M ) + ^F(M c )}0 -e = 0. (31) 
o 2 5 3 
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Performing the Matsubara frequency sums as in Ref.[23, the logarithmic integral which appears in the above 
expression for the effective potential Eq. (|3(JI) divides into two parts. A zero-temperature part Qo(M) which is 
divergent and a finite-temperature part Qp(M) which is finite. For S = M a , Mq , we have 



Q(S) = - / ln(<4 + k 2 + S 2 ) 



Qo(S) + Qp(S) 



d 3 k W£ 



d 3 k 
(2^)3 



ln[l 



(32) 



where tor — \/ (k 2 + S 2 ). For S — M c , we get 



Q{S) = ln{(u> n + im) 2 + k 2 + S 2 ) 



Qo(S) + Qp(S) 



d 3 k ujj: 
(2tt) 3 T 



d 3 k 
(2tt) 3 



[ln(l 



-/3(W+/xr) 



ln(l 



_ e -/ 3 ("fc-A'/)l 



(33) 



Similarly, the second type of integral in Eq. (|29|l is divided into a zero-temperature part Fo(M) and a finite- 
temperature part Fp(M). For S — M a ,M , we have 



F(S) = ( 
Jd 



p to 2 +k 2 + S 2 
F O (S)+F (S) 
d 3 k 1 



(2tt) 3 2u;t 



d 3 k 1 



(2tt) 3 lu^ e^s 



(34) 



And for S = M c , we get 

F(S) 



1 



(w„ + i/i/) 2 + k 2 + S 2 



F (S)+F f3 (S) 
d 3 k I 



(27r) 3 2wr J (2tt) 3 2wr 



d 3 fc 1 



1 



f 



(35) 



The second term vanishes at zero temperature while the first term survives, but the evaluation of the integrals 
in Eqs. I(32|) JSSt IpS t and ((35(1 requires renormalization. Renormalisation in many-body approximation schemes is a 
nontrivial procedure, but in the case of linear sigma model it has been shown by Rischke and Lenaghan [sLllSj. that at 
least in the Hartree approximation the results are not affected qualitatively. We therefore simply drop the divergent 
terms and keep only the finite temperature part of the integrals. A similar procedure has been adopted in a number 
of other investigations using the linear sigma model 

US- 
Using the above definitions of the integrals in Eqs. (|32|) (|33|H|34|I and 1)35(1 . we can rewrite the finite temperature 
effective potential given in Eq. ll.'SOII in a more compact form 

V(4>,M) = ^m 2 4> 2 + ^ 4 -e^ + Q f3 (M a ) + Q p (Mo) + Q,3(M c ) 

-±{F fj (M„)} 2 - ^[F f3 (M )} 2 - ±lF fj (M c )} 2 - ±F p {M a )F p {M ) 



~Fp[M a )F p [M c ) - ^F P (M )F P (M B ). 



(36) 



This final expression contains only terms which are finite. 
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FIG. 3: Solution of the system of gap equations in the case of e / 0. The effective masse of the neutral pion Mo as a function 
of T corresponding to m = OMeV, m = 50MeV, fii = lOOMeV and fn = 138MeV respectively. 
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FIG. 4: Solution of the system of gap equations in the case of e ^ 0. The effective masse of the charged pion M c as a function 
of T corresponding to /u = OMeV, /n = 50MeV, fii = lOOMeV and fn = 138MeV respectively. 
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FIG. 5: Solution of the system of gap equations in the case of e 0. The effective mass of sigma M CT as a function of T 
corresponding to = OMeV, fii = 50MeV, /i/ = lOOMeV and fii = 138MeV respectively. 
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IV. NUMERICAL RESULTS IN THE HARTREE APPROXIMATION 



We have solved the system of Eas. (|28|) and (|31|l using a numerical method based on the Newton-Raphson method 
of solving nonlinear equations. In this way, we are able to determine the effective masses M a , M , M c and the order 
parameter <p as functions of pi and T. We are interested in both negative and positive pi. The critical isospin 
chemical potential at zero temperature is defined by pj — = 138MeV. At T = 0, for \pi\ < m^, the system 
is in the same ground state as that at pi = 0. However while it is favorable to have a nontrivial vacuum which is 
characterized by the appearance of a condensate < 7r+ > or < tt~ > when \pi\ exceeds 0,03,122]] ■ Since in this 
work, we only consider the effects of the isospin chemical potential on the masses of mesons and do not discuss the 
pion condensate 

Ems 

or the pion superfluidity [28j], we constrain our discussion to the values of pi < 138MeV. 

We show in FigOl Fig0]and FigO^the effective masses of the mesons Mo, M c and M a as functions of T corresponding 
to pi = OMeV, pi = 50MeV, pi = lOOMeV and pi = 138MeV respectively Forpi = OMeV, at finite T we reproduce 
the well known results for the effective masses Mo, M c and M CT , i.e., see Refs[l0j. For pi ^ OMeV, the behaviors of 
the meson masses at finite T are very similar to the case of pi = OMeV. At low temperature, the pions appear with 
the observed mass, but their masses increase with temperature while the sigma mass decreases with temperature. At 
high temperature (higher than ~ 300 MeV), due to interactions in the thermal bath, the interactions in the thermal 
bath, all particles (pions and sigma) appear to have the same effective mass. At T = OMeV, the masses of the pions 
and sigma particle remain constant for different pi . 

In FigOll the results for neutral pions at various finite isospin chemical potential are nearly identical to each other. 
This means that the finite isospin chemical potential pi almost can not affect the effective masses of these particle 
duo to the fact that the neutral pions do not contribute to the finite isospin chemical potential pi. However, from 
Eq. (|22c() , for a non- vanishing isospin chemical potential, the masses of the charged pions at the tree level will acquire 
a nontrivial isospin chemical potential dependence for finite values of temperature. In Fig0] as the temperature 
increases, the effective mass of the charged pions at finite temperature changes considerably for different values of 
the isospin chemical potential, especially at the range of T — 50 ~ 300M eV. Moreover, from Fig[3] and FigQ] for a 
certain temperature region, the masses of pions (both charged pion and neutral pion)increase with the increase of the 
isospin chemical potential pi, though the mass of the neutral pion changes very slightly. 

In FigO the results for sigma particle are very difference from the poins. When the temperature is T < 235MeV, 
for a fixed temperature, the mass of sigma decreases considerably with the increase of the isospin chemical potential 
pi, while for temperature T > 235MeV, for fixed temperature, the mass of sigma increases very slightly increases as 
isospin chemical potential pi increases. In this case, the results for the sigma particle mass are nearly identical to the 
ones for the neutral pions. 
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FIG. 6: Evolution of the order parameter <j> as a function of T corresponding to [ii = OMeV, [ii = 50MeV, /ij = lOOMeV and 
fii — 138MeV respectively. 

The presence of the finite isospin chemical potential /i/ modifies the evolution of the order parameter (j>. As it 
is shown in FigEJ at low temperature, for a fixed temperature, the order parameter decreases considerably as the 
isospin chemical potential increases. But at very high temperature, the mass of sigma is very slightly changed with 
the increase of the isospin chemical potential fii. For various pi, when the temperature increases, the order parameter 
decreases, and vanishes smoothly at very high temperatures. In this case, the change is not a phase transition any 
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more and we encounter a smooth crossover from a low temperature phase to a high temperature phase. We conclude 
that the presence of a finite isospin chemical potential when < m, does not change the nature of the chiral phase 
transition, only modifies the effective masses of mesons in the linear sigma model, in accordance with earlier studies 

V. CONCLUSIONS 

In summary, we have discussed the effect of the finite isospin chemical potential to the effective masses of the measons 
at finite temperature in the framework of the 0(4) linear sigma model with explicit chiral symmetry breaking. We 
have provided an alternative mechanism to include the isospin chemical potential in the linear sigma model. In order 
to get some insight into how the the isospin chemical potential could affect the chiral phase transition, we have solved 
the system of gap euqations numerically and found the evolution with temperature of the effective thermal masses of 
mesons at various /ij. 

For neutral pions, the results for different values of finite isospin chemical potential are nearly identical to each other 
due to the fact that the neutral pions do not contribute to the finite isospin chemical potential But for charged 
pions, the situation is very different. At low temperature, the effective masses of charge pions are very different than 
that of neutral pions due to the effects of the finite isospin chemical potential, it appears that they increase more 
dramatically than that of neutral pions as the temperature increasing. But for high temperature, the situation has 
been changed, the effective mass of charged pions are similar to that of neutral pions. On the contray, the results for 
sigma particle are very difference from the poins. For a fixed temperature, when T < 235MeV, the mass of sigma 
decreases considerably with the increasing fj,i, while for T > 235MeV, the mass of sigma increases very slightly as 
Hi increases as well. The evolution of the order parameter (f> for different /i/ shows the effects of the finite isospin 
chemical potential do not change the nature of chiral phase transition, but just modify the effective masses of mesons 
in the linear sigma model. 

Of course, as we have already pointed out, we have not considered the pion condensate 0, |2(| |27| or the pion 
superfluidity psf . Our discussions for are restricted the values of /i/ < 138MeV. Since there will be more complicated 
and richer vacuum structures when \i\ is larger than 138MeV, it is of interest to extend our work to finite temperature 
for /ii > 138MeV and discuss the pion condensate and pion superfluidity. According to the work of Shu and Li, our 
work can be used to study the Bose-Einstein condensation and chiral phase transition with explicit chiral symmetry 
breaking. All these works are in progress. 
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